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Abstract 



We study bound object formation in a background universe accelerated by 
cosmological constant and by quintessence. If the acceleration lasts forever, 



> 

en 

Tij- . due to the existence of event horizon, one would have naively expected the 

^— ^ . universe to approach a state of cold death. However, we find that many 



local regions in the universe can in fact be protected by their own gravity to 
form mini-universes, provided that their present matter densities exceed some 
critical value. In the case with cosmological constant (ACDM cosmology) 
the condition of forming mini-universe is that the ratio of present density 



^ . parameters il^/i^A ought to be larger than a critical value 3.63. Such mini- 

^ universes typically weigh less than 2 x 10^^ solar masses, with the lighter ones 

having tight and compact configurations. In the case with quintessence, the 

final ratio Qi/i^q of a mini-universe is found to be always larger than w'^ — Wg, 

where Wg is the equation of state parameter. 
PACS numbers: 98.80.-k, 98.80.Es 
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I. INTRODUCTION 

Recently direct evidences from the studies of Hubble diagram for Type la supernovae |I| 
indicate that our universe is expanding with an increasing rate, i.e., accelerating expansion. 
This result compels one to seriously consider a dominant component of exotic dark energy, 
the quintessence, in the cosmic energy balance. Accelerating expansion implies that the 
deceleration parameter q = {Qm + (3wq + l)Qq)/2 to be negative, where Wg is given by the 
cosmic equation of state, Pg = WgPg, of the quintessence Q. As long as Wg < —1/3, the 
corresponding dark energy provides a negative contribution to q, and the relevant range of 
Wg for acceleration lies in between —1/3 and —1. When Wg assumes the extreme value —1, 
constrained by the weak energy condition 0, the dark energy is reduced to the one arising 
from the cosmological constant A. 

The acceleration may stop in the future, i.e., a transient phenomenon with a time- 
dependent Wg, or may last forever, depending on the nature of the dark energy 0. A positive 
cosmological constant, consistent with the supernovae data aX z = 1.7 [Q, leads to a forever 
acceleration. Quintessence scenario with a constant Wg, not ruled out by data, can also lead 
to forever acceleration. There are profound implications for a forever accelerating universe 
PJ^. The universe will exhibit an event horizon; that is, there exist regions of the universe 
inaccessible to light probes. It has been argued that such a universe presents a challenge 
for string theories or any of its alternatives f^. Many other implications of cosmological 
constant and quintessence and specific models have been studied in the literature [0-§. 

Naively speaking, in a forever accelerating universe, two points at different spatial lo- 
cations will, as measured at any one of the two points, eventually be pulled apart at the 
speed of light to approach the future horizon. Hence existing structures will be eternally 
frozen and the universe will approach a state of cold death. Is this an unavoidable conse- 
quence everywhere in the universe, or may some mechanisms be in operation that salvage 
the future? The present work confronts this issue and shows that the self-gravity of matter 
can protect a local universe from ending up with cold death, provided that the present local 



matter density is sufficiently high. Contrary to the above naive expectation, a large number 
of mini-universes, including our own Local Group, will become self-bound and survive the 
destruction of cosmic repulsive force due to the quintessence. They can naturally arise de- 
spite that these mini- universes will still be falling into each other's horizon at a sufficiently 
late time and become isolated island universes. 

At present, a variety of dark energy candidates can accommodate the experimental data. 
If the dark energy arises from a positive cosmological constant, the AC DM cosmology, 
data from recent cosmic microwave background (CMB) radiation |]lOl further constrain the 



present energy density parameters to be Q^ = 0.35 and 1)^ = 0.65. While for quintessence, 
depending on the value Wg, the present matter density fl^ and quintessence density i7° can 



be different 11 



For our purpose of addressing the future of local bound objects in an accelerating uni- 
verse, the cosmological constant presents the worst scenario since the cosmological constant 
gives rise to the strongest repulsive force against formation of bound objects. If objects are 
not torn up by the cosmological constant repulsive force, they will survive in other scenarios. 

The paper is arranged as the following. In Sec. (2), we ffist consider the ACDM cosmology 
in details, by assuming that the present universe is spatially fiat with 65% energy given 
by the cosmological constant {Q% = 0.65) and another 35% by the matter (fi[^ = 0.35). 
In Sec. (3), we then consider the case where the cosmological constant is replaced by the 
quintessence with a constant Wg. In both cases the mini-universe can form by self-gravity, 
but the properties of critically bound mini- universes are very different. Our conclusion is 
given in Sec. (4). 

II. BOUND OBJECT FORMATION WITH COSMOLOGICAL CONSTANT 

We shall ffist adopt a toy model to address the questions of whether bound objects may 
still form after the background expansion has started to accelerate and how these objects 
become virialized. We then study a more realistic model where the observation/simulation- 



motivated features of virialized objects are taken into account. 

A. Self-Binding by Local Gravity 

The toy evolutionary model for bound object formation assumes the universe to consist 
of the background and an isolated cold dark matter fluctuation that grew from a small- 
amplitude since the early epoch. The localized perturbation possesses spherical symmetry 
and consists of two concentric solid spheres of different densities. The inner sphere has a 
radius Tj with a uniform dark-matter density pi^ which is greater than the background matter 
density Pm- The outer sphere has a radius Tout and a uniform compensating under- density, 
such that the averaged matter density within Tout equals pm- Any particle outside Tout feels 
no extra gravity resulting from the over-density and expands as it would have been in a 
homogeneous universe. The inner sphere can be regarded as a closed Friedmann universe, 
which has a Hubble parameter. Hi, given by 



rA^ SttG Ki 



H^=[i) =^(p^ + Pa)-3, (1) 



where Pa = (l/87rG')A is a constant and /tj is the curvature of the mini-universe. 
A test particle at Tj feels an effective potential V{ri), 

,,, , GMi AttG 2 

Viu) = '- - ^-pAvl (2) 

Conservation of matter yields a constant mass Mj within Tj with Mi = 47rpi(rj)rf/3. The 
potential V{r) has a maximum at Vmax determined by V'{rmax) = with V{rmax) = 
—AnGpAr'^ax a-nd Vmax = {p'i/'2pA)Y^^i"Q, where the density at present has been denoted 
as p° = Pi{ro) with tq being the present radius of the over-dense sphere. When the test 
particle reaches Vmax with a zero velocity, it will be marginally bound. The density pi cor- 
responding to this situation is called the critical binding density, denoted as p^^ at present, 
and any local region with density larger than p°^ will eventually turn around and collapse. 
The curvature of the mini-universe can be conveniently obtained at the turn-around. 



after which gravitational collapse ensues, by setting r,^ = and equating —Ki/2 to the 
potential at the turn-around. That is, 



rCi 



3 PA rf„ 
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where r^ is the turn-around radius. In particular for a critically bound mini-universe, we 
have rta = Tmax and Kic = Ki{rmax) = STrGpAr^^^. 

The time span needed for the mini-universe to evolve to the present epoch is given by 
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(4) 



VG^^PA Jo v'(p?/Pa) +y^- (3«:,/87rG'pAr2)y2/3' 
When ti is known, Eq.(4) can be inverted to solve for the present local density p°. Thus one 
needs to obtain information about tj. To this end we note that the mini- universe should have 
started to grow from a small amplitude perturbation at the beginning of matter domination 
around 2; = 3 x 10'^. Therefore to a good approximation, the evolution time ti is equal to 
the background evolution time Ib- The time ts for the background Friedmann universe to 
reach the present ratio Vt^/Vt\ is 



/■«o da 1 /■! dy 1 . , -i 

ts = -^ = 1^ ^ I I = 1^ ^ smh ' 

Jo a 



"- (5) 
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Setting ti = tB, one can determine the present density ratio p°/pA for a given Kj. 

We consider three cases of interest for the discussions of the density ratio Pi/pA'- i) the 
critical mini-universe, where the test particle reaches Vmax with zero velocity and is therefore 
marginally bound; ii) the case where the over-density is larger than that of case i) and the 
test particle has a zero velocity and turns around at present; iii) the over-density is even 
larger so that the turn around occurred at the time when the background matter density 
p'^ equaled pa, a particular epoch prior to which the matter in- fall had been vigorous. 

For case i), we have /tj = /tjc = S^rGpA^maa;- Given the background density parameters 
il^ = 0.35 and ^2^ = 0.65, Eq.(4) along with Eq.(5) yield the present critical density ratio 
P?c/Pa = 3.63 {Qic{ro) = 2.36) and ro = rmax/l-'22. 



We note that over-dense regions with the present density parameter larger than the 
critical value r2jc(ro) are abundant in the universe. In fact nearly all future bound objects 
will never reach Vmax with a turn-around radius rta smaller than Vmax- Inserting k, of Eq.(3) 
into Eq.(^), one obtains their present densities pi{ro) for a given ratio rta/ro. For case ii) 
the mini-universe that is presently turning around, we have rta = fo, and it is found that 
pV Pk = 5.90 with fij^o = 3.83. Any over-dense region with a density greater than this value 
has already undergone collapse. For case iii), it is found that pT I Pk = 8.55 with Qi^eq = 5.52. 

Although the toy evolutionary model so far discussed may seem somewhat simplistic, the 
model nevertheless illustrates the essential physics ||12| and supports our contention. From 
the above discussions, we clearly see that local bound objects can form against the cosmic 
repulsive force, as long as the local density parameter fii(ro) at present is greater than 2.36. 
This value is at least a factor of several less than the mean density parameter of the Local 
Group, which includes the Milky Way, the Andromeda galaxy and more than a dozen of 
smaller galaxies such as the Magellanic clouds. It suggests that the Local Group has defied 
the cosmic repulsive force and gravitational collapse is underway. 

B. Virialization 

Almost all bound objects, except for a set of measure-zero, critically bound ones will 
undergo, or have already undergone, gravitational collapse. It is instructive to further 
examine the issue of virialization after an object collapses. The virial density found 
here helps us address deeper questions concerning the internal structures of bound ob- 
jects to be addressed below. The virialized kinetic energy of a test particle is given by, 

< T >= — < F ■ f > /2. For the simplified model under consideration, < T > is given by 

< T >= {GMi/2ry) — (AnGp/^rl/S), where r^ is the virialized radius. This leads to 

-^=<T> + <V>= -(^ + ^-pAr^). (6) 

From both energy and mass conservation together with the expression of ki at the turn- 



around, eq.(H), with the subscript "0" replaced by "ta", one finds [|T^] 
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Pi^^2 - —) + 2{1 - 2^) = 0. (7) 

PA r^ rfa 

Regions with a space curvature k, infinitesimally greater than Kjc of case i) will finally 
collapse in the infinite future, and we shall also regard this situation as case i) in a broad 
sense. It is straightforward to obtain that r^/rta = 1/2.73 for case i), since pi.tal Pk = 2 and 
the corresponding r^/rQ is 1/2.24. The ratio of the virialized density to the A energy density 
is then found to be Pv/pa = 41. For cases ii), we find r^/rQ = 1/2.20 and Pv/pa = 62. For 
case iii), r^/rta = 1/2.13 and p^/pA = 83. 

Case (iii) is interesting from the consideration that the matter in-fall was considerably 
vigorous prior to the turn-around but became suppressed afterward. Although its collapse 
started aX z = 0.23 when Vim = ^a, as the completion of collapse takes a time span twice the 
turn-around time, this object will be virialized in the future when the background universe 
reaches ^a = 0.89. By then the background scaling factor a{t) becomes a factor 1.63 times 
the present value ag. 

After the bound object collapses and forms a virialized object, what does the bound 
object need to adjust so as to guard itself against the stripping repulsive force? We shall 
address this question by a model of somewhat more sophistication. It has been empirically 
found by cosmological A^-body simulations, which compute at best up to the present epoch, 
that cold dark matter particles tend to cluster into individual bound objects with a universal 
density profile. This density profile has been fitted by the empirical formula: 

P = ( 1 ^2 ' (8) 

ryVs + ry 

where A and r^ characterize the mass and inner scale length of the bound object |]T1 . 



This density profile has a r~^ cusp and a r~^ envelope extended out to the virial radius 
r„. It is universal in that the parameterized functional form is independent of cosmological 
parameters and epochs, although the coefficients A and r^ are. Outside the virial radius r^, 
the density profile must be truncated so that the interior matter can be blended into the 
background matter. Typically the concentration parameter t^/ts exceeds 5 for the present- 



day virialized objects, and the r~^ envelope in the density profile is noticeable. Such an 
universal profile has been tested against observations favorably | |15|| . 

For our purposes, the universal profile can be extrapolated to describe the configurations 
of future bound objects. We will be concerned about how the concentration parameter r^/vs 
and the characteristic density A vary with the bound-object formation time. Unfortunately, 
the virial radius r^ is undefined in Eq.(||), therefore prohibiting us from using this density 
profile directly. One hence needs to seek another method to address the evolution of the 
universal profile. 

The following strategy may be adopted for a quantitative evaluation of the evolution 
of the profile. The steady-state collisionless dark matter is described by the phase-space 
distribution function /(Ci(x, p), C2(x, p), ...), where Cj's are the constants of motion. For 
convenience, we assume the momentum p to be isotropic and the mass density spherically 
symmetric. It follows /(Cj) = f{E), a function of only particle energy E. Set f{E) = 
for unbound particles that have E > 0. Upon integrating out the momentum from f{E){= 
f{p'^/2m + 0)), we obtain the mass density p(0), a function of the potential of all forces. 
The universal profile Eq.(|^) is obtained from the Poisson equation for an appropriate p{4>). 
The functional form of p(0) can be approximated by the following considerations. The inner 
singular density profile requires p(0) ~ (0 — 0mm) near the potential minimum 4>min at r = 0. 
On the other hand, the outer r~^ behavior requires p(0) ~ —0'^ at the edge of the potential 
well, — > 0, where the bound particles are about to escape. Thus these desired behaviors 
can be captured by 

P^'^^ = "exp[/5(0-0^,„)]-l- ^^^ 

Normalizing all quantities in dimensionless forms, we arrive at the Poisson equation, 

x"^ dx dx exp[6(?/^ + 1)] — l' 



where ip = 0/|0mm|(< 0), x = r/L with the characteristic length L = y |0mm|/87rG'pA, 
b = P\4>min\ and g = Q;|0mmP/2pA- The gravitational potential 0c, on the left-hand side of 
the Poisson equation has been replaced by + {4:iTG/3)p\r'^. 
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The solutions to eq. ([T0|) require numerical integration. Different values of g and b yield 
different A and r^ of eq.(§). This is a nonlinear eigenvalue problem, and there exists a locus 
in the eigenvalue {g, h) space, along which the bound objects satisfy the desired boundary 
condition that p vanishes at the potential maximum. Figure (1) presents the parameter 
6 as a function of g for these bound objects. It is instructive to seek a measure for the 
the inner r~^ core region. Also plotted in Fig.(l) is the dimensionless half-mass-weighted 
radius h = {xD)i/2 of these virialized objects, where D = p/2pA- The quantity {xD)i/2 = 
3 Jq^^^ x^Ddx/xy2, and xi/2 is defined to satisfy Jq^ x^Ddx = (1/2) Jq""""' x^Ddx with 
D{xmax) = 0. Note that the parameter g varies by a factor of 10^ and b by several hundred, 
and yet the quantity h does not change in most parameter space until near the high end of 
{g,b), where it decreases only by 50%. This result shows that the r~^ core contains more 
than half of the mass for most solutions, except for those with g > 5 x 10^. 

Plotted in Fig. (2) are the density and potential profiles of three solutions with {g,b) = 
(12, 0.0063), (5x 10^ 2), and (8x 10^ 3.2). These solutions, denoted as (a), (6) and (c) respec- 
tively, have progressively stronger cores of smaller sizes as both g and b increase. Solution 
(a), having an extended r~^ core, gives no detectable outer envelope in the density profile, 
whereas solution (c), with a strong and compact core, contains a distinct r~^ envelope. 
These solutions are expected to have progressively larger ratios of Pv/pA, corresponding to 
virialized objects formed at progressively earlier epochs. 

As virialized objects formed at different epochs appear differently, a future observer ought 
to be able to date these objects by morphology. To determine the formation epochs, we first 
start from the extended-core solution (a), which is at the small end of the {g,b) spectrum 
shown in Fig.(l). It has already been well within a self-similar scaling regime beginning 
around {g, b) ~ (196, 0.1), where the denominator of p(0) in Eq.(^ can be expanded to the 
leading order of small /5, yielding g/b -^ 1953. Different solutions in the scaling regime have 
the same profiles and hence solution (a) can be identified to be the virialized configuration 
of the last collapsed objects. Making use of the result derived earlier that Pv/ Pk = 41 within 
the virial radius x„ for the last bound object, we determine the mass fraction, equal to 0.972, 



within the virial radius. 

Assuming that the mass fractions, 0.972, within the virial radii of all solutions remain 
the same, we find that pv/ Pk = 65 and 85 for solutions (b) and (c), respectively. According 
to the virialization model analyzed earlier, solutions (b) and (c) correspond to the virialized 
objects with turn-around times at 2; = and z = 0.23, or cases ii) and iii) considered earlier, 
and they are to complete the collapse in the future when the scaling factors a/ao become 
2.25 and 1.63, respectively. 

We note that the redshift evolution of morphology can be rather rapid for these future 
virialized objects formed in between a/ao ~ 1-5 and 4.5; by contrast, the collapsed objects 
up to the present do not show much difference in morphology, all having a compact core 
and a r~^ envelope. 

Common to all virialized objects shown in Fig. (2) is that matter is mostly confined at 
a considerable distance away from the vacuum radius Xmax- Moreover, as the potential i/) 
vanishes quadratically (x — XmaxY near the potential maximum, the density should vanish 
as (x — Xmax)^, accordiug to eq.(|^). Using the force balance of pressure against gravity, it 
means that the local temperature vanishes as (x — Xmax)"^ at the edge. That is, the dark 
matter is cold at the outer edge of the mini- universe. Moreover, conservation of the phase- 
space volume for collisionless particles demands that Ar^Av^ remains a constant, and when 
Av^ — !► we have Ar —>■ 00. That is, the mini- universe contains a thick outer layer of cold 
dark matter particles. 

How do these cold particles result from? Some fraction of collapsed dark-matter particles 
must have been pulled away by the cosmic repulsive force to the background during virial- 
ization. Thus, evaporation should serve as an important cooling mechanism to produce a 
thick outer layer of cold particles. The mass loss due to evaporation may be estimated. We 
note from Fig. (2) that the cool matter outside the virial radius r^ occupies about 90% of the 
volume enclosed by Xmax, but it amounts to only 2.8% of the total mass. The evaporated 
material can not be much more than the cool matter located in the thick outer layer. There- 
fore the small mass loss should not create any problem for the above quantitative analyses 
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that assume mass conservation during virialization. 

Finally, the typical masses of bound objects formed at various epochs are also of rele- 
vance toward understanding the future world. As the linear power spectrum P{k) is well 



understood |T6| , p^ , this issue can be properly addressed as well. With a given conserved k,, 
the radius r, can be traced back to a sufficiently high z through Eqs.(^) and (|^), thereby 
fixing the amplitude of linear over-density by the relation Sp{z) = 3K,i/8TTGrf{z). On the 
other hand, Sp{z) is also given by g{z) Jq A7ck^P(k)dk, where g(z) is the growth factor 

of linear perturbations [|l^. Identification of the two (5p's allows ri{z) to be determined, 
and the enclosed mass is simply M = A'n-p{z)rf{z)/3, a redshift-independent quantity, as it 
should. The typical mass of the last bound objects (solution (a) of Fig. (2)) is found to be 
2 X 10^^ solar masses, the mass of a poor galaxy cluster, and those for solution (b) and (c) 
are 10^^ and 4.5 x 10^^ solar masses, respectively (assuming Hq = 65 km/sec/Mpc). 

III. BOUND OBJECT FORMATION WITH QUINTESSENCE 

In this section we study the condition for bound object formation under the repulsive 
force of the quintessence with a constant Wg. At present experimental data allow for the 
possibility that the accelerating cosmic expansion be due to quintessence with a relatively 



wide range of constant Wg{> —1) and energy densities |TT|]. 

In order to have an accelerating universe at present, the equation-of-state parameter Wq 
should satisfy 

Unfortunately, the value of Q^/fi'^ is still rather uncertain, and its plausible value deduced 
from data correlates with Wg, with a larger Q^/Q^ for a smaller Wg in a flat universe. For 
example, when fi^ = 35% and fi^ = 65% to make up the energy budget for a flat universe, 
eq.(|TT]) implies Wg < —0.5, but the existing data pushes the plausible Wg to be near —1. 
Nevertheless, to have an idea as to how the results vary with Wq, we shall consider three 
representative values, Wq = —0.5, —0.75 and —1, in our later discussions. 

11 



The quintessence energy density varies with the scahng factor as Pq = p°(a/ao) 
The reference background cosmological time ts in this case is given by 



-3il+Wq) 



tB = J- C '^^ =. (12) 

For the study of bound mini- universes, the situation becomes more comphcated than 
that with a cosmological constant because the space curvature Kj within the collapse region 
is time-dependent, and eq.(4) is no longer valid. (The analysis presented in Ref. [|T^], which 
assumes a constant Kj, is hence invalid.) This problem should be analyzed using the momen- 
tum component of the Einstein equation, which describes the force balance of inertia with 
gravity and quintessence force. It can be understood as follows. Since the self-gravity of 
quintessence is negligible, the repulsive force of quintessence at a radius r varies with time as 
{a^l a{t)f'^^^'^''' . This is a time- dependent force, thereby leading to a time-dependent space 
curvature. 

The momentum component of the Einstein equation is a second-order differential equa- 
tion: 

f AtxG , 



r 3 

Defining x = {a/ao), y = (r/ro), one has 



[p, + {l + 3wq)pq]. (13) 



i = \/^p°v'(fiyfi°)(i/^) +^"(^+=^"'^^ 

/I (^ r)0 

y = -^PUTToy'" + (1 + 3^,)x-^(^+-^)]. (14) 



3 ^^"^fio^ 
Combining these two equations one arrives at 



2y"xil + x-'^""') = y'(l + (1 + 3wq)x-^'"') 

— 2 — 



+ {l + 3wq)!^x-'-^), (15) 

y^ X 



where x = a;((]^/f]°)i/3«;,^ y ^ ^((]0^/nO)i/3«,,^^o^/^0) ^^^ g ^ ^^/^^_ 

We may solve eq.(15) to obtain the conditions for bound object formation at any back- 
ground cosmic time ts, c.f., eq. ([T2|) . However, due to the wide range of possible parameters 
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for the quintessence cosmology, we shall instead concentrate on the critically bound mini- 
universe. The asymptotic state of the critical mini-universe can be found by letting both 
y —>■ oo and a; ^ cxo in the infinite future. The existence of such a state also requires that 
the background expansion be accelerating, i.e., Wq < —1/3. It then follows that eq.(15) has 
a scaling solution, y = {Qiao/^qooY^^x^^^'^''\ where 



^2 



Wq - Wq, (16) 



by matching the coefficients in eq.(15). The subscript "oo" refers to evaluation of a quantity 
at a reference epoch in the infinite future. 

The space curvature Kj for this mini-universe defined as minus of twice the sum of kinetic 
energy and potential energy (of the force used in eq.(p!4D) is given, in the asymptotic limit, 

by 

n^ = -(— )-(i+-^)(l + 3Wq)AnGpq^rl, (17) 

where eq.(^) has been used. This confirms the statement made earlier that with 
quintessence the curvature is time dependent. 

The space curvature decreases with time as a~^^^'^''^ but its magnitude cannot be deter- 
mined from the scaling solution; it can be fixed only by integrating eq.(15) to obtain the 
entire solution. In this asymptotic state, the absolute values of kinetic, potential and total 
energies all have the same time dependence with comparable magnitudes. For Wq = —1, we 
recover the cosmological constant case with pioo/Pgoo = 2 and Kj given in eq.(|D. The above 
critical mini-universe is very different from that of the case with a cosmological constant, in 
that its size grows indefinitely since r grows as a^^^'^''\ 

One may further extrapolate the ratio Pioo/Pgoo back to the present by integrating eq.(15) 
numerically. The results are shown in Fig. 3. In Fig. 3 we show log{Qi/Qq) as a function 
of \og{Qq/Qm) for Wq = —0.5, —0.75 and —1. From Fig. 3, one can read off the densities of 
the critically bound mini-universe on the vertical axis for various background densities on 
the horizontal axis. (The horizontal axis can be regarded as another way for expressing the 
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background cosmic time, c.f., eq.(|T^).) Regions whose present local densities exceed by a 
finite amount above a given curve shown in Fig. 3 will turn around within a finite time and 
become virialized in a manner similar to that discussed in Sec. (2), though the details differ 

However, it is interesting to note that there are differences for the virialized objects 
between the general quintessence cosmology and the ACDM cosmology. For the case with 
quintessence, a local region with an energy infinitesimally smaller than the critically-bound 
condition (local matter density larger than the critical density by an infinitesimal amount), 
the region will undergo collapse in the infinite future. The bound object has a nearly zero 
binding energy, as can be seen from eq. (p!?!), and hence after virialization its virial radius 
Vy discussed in Sec. (2) will be infinitely large, a great contrast to the Wg = —1 case where 
the just bound objects have finite sizes. 

As has been mentioned earlier in this section, the precise value of fl^/Q^ is yet to 
be determined and the plausible value of Q'^/Q^ in fact anti-correlates with the value of 
Wg. Though for Wg = —1 the most plausible fij^/fi^ = 0.35/0.65, for Wg = —0.75 and 
—0.5 the most plausible fi^/fi'^ decreases to approximately 0.1/0.9 and (^)0/l, respectively 
|rT|. Precise knowledge for the content of energy forms, which can be acquired from more 



accurate calibration of Type la supernovae and measurements CMB radiation as well as 
other observations [0,0, are crucial for making good use of Fig.3. 

IV. DISCUSSIONS AND CONCLUSIONS 

In this work, we have quantitatively shown that mini-universes can form against the re- 
pulsive force of cosmological constant or quintessence with a constant equation-of-state pa- 
rameter Wg, as long as they have sufficiently high densities at present. These mini-universes, 
when virialized, have distinctly different configurations depending on the formation time. 
Bound objects form earlier appear more compact and are confined by deeper gravitational 
potentials. This issue is addressed in details for the ACDM cosmology, where the mini- 
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universes are found typically no heavier than a poor galaxy cluster. 

As discussed in Sec. (2), the actual size of mini-universe varies, depending on the pri- 
mordial density power spectrum. It is instructive to examine mini-universes of very 
large scale. We note that a marginally bound mini-universe has a gravitational radius 
r^ = c^r^ax/'^i — c^/(3^A-f^o) i^ ^^^ '^^^^ ^^ A-cosmology. If a marginally bound mini- 
universe has a physical size larger than Vg, no light emitted within can escape from it, and 
hence behaves like a black hole when observed from the exterior. Such a mini- universe is 
absolutely stable against the background pulling force. We therefore obtain the minimal size 
i^min of an absolutely stable and marginally bound mini-universe to be c/\^3Q\Ho, which is 
comparable to the Hubble radius of the background universe. 

For mini-universes with larger local densities, the required radius is smaller. These mini- 
universes will eventually undergo catastrophic collapse, and can asymptotically be described 
by a static space-time metric - the Schwartzschild-de Sitter metric: 

d.- (1 - ?^ - Ir'^H - ( ,_,c^;:_^,.^3 + r^d'O + r' s.u^ M',). (18) 

It has two horizons: an inner horizon at a modified Schwartzschild radius and an outer 
horizon pertaining to the accelerating cosmic expansion. The region in between the two 
horizons is unstable, in that an observer located in between the two horizons must eventually 
fall onto either of them. 

Our local universe cannot be within such a closed universe, as the angular scale of the 
first peak in the CMB anisotropy power spectrum indicates the local universe within the 
Hubble radius to be spatially flat [|1^. This, however, does not preclude the possibility 



that regions outside the present Hubble horizon may have sufficiently large super-horizon 
over-densities to become a closed universe. The required present matter over-density Q'^/Q^ 
in such universes is as high as 6.75. These regions are very rare if the primordial density 
fluctuation is indeed scale- invariant, created during the inflation in early universe, with a 
power spectrum linearly proportional to the wavenumber k. Such local closed universes 
should be more abundant in the past, and some of them may have been within our present 
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horizon, perhaps in the form of super-massive black holes. 

In the case of quintessence, the situation is more complicated because the mini-universe 
has a time-dependent space curvature. The curvature actually decreases with time, due to 
the fact that matter and quintessence as a whole can be regarded as a non-ideal fluid and 
the local matter over-density gives rise to entropy fluctuations. Space curvature fluctuations 
pertinent to entropy fluctuations always decrease monotonically in time, regardless whether 
the fluctuations are of sub-horizon size or super-horizon size. However, adiabatic super- 
horizon fluctuations can also give rise to curvature fluctuations with no time dependence, 
for which quintessence fluctuations can grow and collapse as well. They are more relevant 
for the discussions of the local closed universe in the context of quintessence cosmology. 
Unfortunately, our analysis presented in Sec. (3) is confined to entropy fluctuations, where 
the quintessence fluctuation is negligible. Nevertheless, the dynamics of such a local closed 
universe with a constant space curvature is relatively simple and like that of a homoge- 
neous universe. Thus the discussions given in the preceding paragraph for A-cosmology also 
holds for quintessence-cosmology. The super-massive black hole so formed should contain 
a substantial amount of dense quintessence field within the Schwartzschild radius, and its 
asymptotic space-time metric depends on the yet unknown nature of the quintessence field. 

Return to the mini- universes arising from entropy fluctuations and of sub-horizon size 
analyzed in Sec. (3). As long as the matter density is larger than the critical density given in 
Fig. 3 for the corresponding equation of state parameter, the local region will be dominated 
by the matter gravity and collapses to become a mini- universe. These mini-universes will 



become virialized in a similar manner as those of the ACDM cosmology [^, due to that 
virialization depends mainly on the balance between gravity and pressure, and to a lesser 
degree on the large-scale cosmic pulling force. Despite the similarity, the details can be 
different. For example, conservation of total energy within the bound object, which is 
essential for applying the virial condition, no longer holds in quintessence-cosmology, and 
hence the resulting virial over-density within the object deviates from that of the ACDM 
cosmology. Also the marginally bound mini-universe has an infinitely large virial radius. 
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At any rate, a single tightly bound mini-universe will, in the future, emerge from the 
Local Group, whose appearance resembles solution (c) of Fig. (2) in the case of ACDM 
cosmology. Interior to the mini-universe, local physics in the future will remain the same 
as what we experience presently. However, the mini-universes will become isolated island 
universes in the future where no communication among them is possible, and each is a 
lone world. When observed within a given mini- universe, all other mini-universes originally 
within the horizon will eventually fall onto the horizon and become frozen. Lights emitted 
from the horizon are frozen as well, making other island-universes dark and invisible. 

The detailed properties of the mini-universes, depend on whether the accelerating ex- 
pansion is due to the cosmological constant or quintessence. New experimental data from 
various observations in the coming decades are expected to reveal more information. By 
then we will have a better understanding of the properties of mini-universes. If the acceler- 
ation of our universe at present is indeed due to cosmological constant or quintessence with 
constant Wq, our universe will enter a new era, the era where all island universes are falling 
onto each other's horizon and appear to fade away. It is, nevertheless, comforting to know 
that the Milky Way galaxy is in one of these island universes, within each of which most 
physics remains the same as is today. Review on the details of how an island universe will 
eventually end itself with physics known today can be found in Ref. [19]. 
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FIGURES 




FIG. 1. Nonlinear eigenvalues {g, h) for Eg. (p^) . When b < 0.1, the solutions are in a similarity 
scaling regime. Also shown is the half-mass weighted radius, h. The constancy of h implies that 
the r~^ core contains more than half of the total mass. 



Density and Potential Profiles 




FIG. 2. Dimensionless densities D (in log scale) and potentials "0 for solutions (a), (6) and (c) 
ofEq.(|lO|). 
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FIG. 3. The barely bound mini-universe critical density as a function of time. The vertical 
axis is log(r2j/f]q) and the horizontal axis is \og{0,q/i^m)- 
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